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ABSTRACT 

Unless properly corrected for, the existence of intergalactic dust will introduce a red- 
shift dependent magnitude offset to standard candle sources. This would lead to over- 
estimated luminosity distances compared to a dust-free universe and bias the cos- 
mological parameter estimation as derived from, e.g., Type la supernovae observa- 
tions. In this paper, we model the optical extinction and X-ray scattering properties 
of intergalactic dust grains to constrain the intergalactic opacity using a combined 
analysis of observed quasar colours and the soft X-ray background. Quasar colours 
effectively constrain the amount of intergalactic dust grains smaller than ~ 0.2 /im, to 
the point where we expect the corresponding systematic error in the Type la supernova 
magnitude- rcdshift relation to be sub-dominant. Soft X-ray background observations 
are helpful in improving the constraints on very large dust grains for which the amount 
of optical reddening is very small and therefore is more difficult to correct for. Our 
current upper limit corresponds to ~ 0.25 magnitude dimming at optical wavelengths 
for a source at redshift z = 1, which is too small to alleviate the need for dark energy 
but large in terms of relative error. However, we expect it to be possible to lower this 
bound considerably with an improved understanding of the possible sources of the 
X-ray background, in combination with observations of compact X-ray sources such 
as Active Galactic Nuclei. 

Key words: cosmology: theory, scattering, (ISM:) dust, extinction, (galaxies:) inter- 
galactic medium, quasars. 



1 INTRODUCTION 

The use of Type la supernovae (SNe la) as standardized can- 
dles to probe the redshift-distance relation remains essen- 
tial for establishing and exploring the dark energy universe. 
The fact that complementary cosmological probes agree on 
the concordance cosmological model with current acceler- 
ated cosmological expansion makes it unlikely that the ob- 
served dimming of SNe la are solely due to dust extinction. 
Nonetheless, it is still possible that dust extinction could 
bias cosmological parameter estimation. Also, when trying 
to explore specific dark energy properties, it becomes crucial 
to be able to constrain the effect of dust o n SN la observa- 
tions l|Menard et al.ll2010J : ICorasanitil [20071 ). 

The presence of dust in the intergalactic medium (IGM) 
has been the subject of numerous studies. Based on esti- 
mates of the stellar density and metallicity as a function of 
redshift, several authors have inferred the existence of sig- 
nificant amounts of dust in the IGM with density fWst ~ 
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10~ 6 - 1CT 5 dLoeb fc Haimanlll997l : llnoue fc Kamavall2004l : 
lFukugitall201ll ). 

Dust grains scatter and absorb photons with an energy 
dependent cross-section and typical dust extinction is corre- 
lated with reddening of the incoming light. The amount of 
reddening is usually quantified by the total-to-selective ex- 
tinction ratio Rv = Av / E(B — V) — Av/(Ab — Av), where 
E(B — V) is the colour excess and Ab a nd Ay are the e xtinc - 
tion in the B and V-band respectively. iMenard et alj (|20ld ) 
report a statistical detection of dust reddening of quasars 
(QSOs) out to large distances (a few Mpc) around galaxies 
at z = 0.3. The observed reddening implies a slope of the 
extinction curve, Rv = 3.9 ± 2.6, which is consistent with 
that of Milky Way dust (Rv = 3.1), albeit with large uncer- 
tainties. Extrapolating this result to higher redshifts yield 
a lower limit estimate of the extinction in the restframe B- 
band of Ab(z = 1) Jj0.03 magnitudes. 

lAguirrd (|l999l ) and iBianchi fc Ferraral |2005l ) find that 
astrophysical processes which transfer dust into the inter- 
galactic medium would preferentially destroy small grains, 
leaving only grains larger than a ~ 0.1/un, implying less 
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reddening and higher values for Rv (often labeled "grey" 
dust). The absence of detectable systematic reddening of 
SNe la with increasing redshift implies that any intergalac- 
tic dust extinction must be quite grey at optical wavelengths, 
suggesting the grains must be large. Since the dust correc- 
tion for SNe la depends on the reddening, large grain s will 
more difficult to cor r ect fo r. [Mortscll & Goobarl 120031 ) and 
lOstman fc Mortselll (|2005l ) simulated the reddening by in- 
tergalactic dust based on the different parametrizations of 
mean extinction laws of Milky Way-like dust, with < Rv < 
12. They used observations of QSO colours and template 
spectra to put an upper limit on the dimming in the rest- 
frame B-band by intergalactic dust of a source at redshift 
z = 1 of A B (z = 1) ^0.02 for R v ~ 3 and A B (z = 1) <;0.1 
for R v ~ 10. 

Apart from the reddening of cosmological sources, the 
presence of dust grains can also be inferred from the ab- 
sorption of UV/optical photons which are re-emitted in the 
far-infrared. lAguirre fe Hainia n (2000) calculated the con- 
tribution from intergalactic dust to the cosmic far-infrared 
background, assuming it was heated by intergalactic radia- 
tion. They found that dust densities of fidust ~ a few x 10 -5 , 
necessary to account for the dimming of SNe la, would pro- 
duce most of the far-infrared background. However, obser- 
vational data leaves little room for any such diffuse emission 
component, since discrete sources detected by the SCUBA 
survey account for alm ost all of the background at 850 /im 
l|Hauser fe Dwekll200ll ). 

Intergalactic dust would also scatter X-rays from point 
sources into extended, diffuse X-ray halos. If the size of 
the halo is large enough, the scattered radiation will ef- 
fectively be part of th e unresolved soft X-ray background. 
iDiikstra fc Loebl (|2009h argue that dust scattered X-ray ha- 
los around Active Galactic Nuclei (AGN) can maximally 
account for a fraction /halo ~ 5 — 15% of the total mea- 
sured Soft X-ray Background (SXB). This allows them to 
place an upper limit on the optical/near-infrared extinction 
Am{z = 1,A = 8269A) <;0.15(/ ha io/10%). The absence of 
an X-ray halo around a sing l e z = 4.3 QSO observed with 
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to place upper limits on the dust density parameter 
fWst S^IO -6 — 10 -5 assuming a constant comoving number 
density of dust grains of size a = 1 /im or with a power law 
distribution of grain sizes 0.1 ^ a ^ 1 ^m. 

In this paper we will combine a QSO colour analysis 
with a SXB analysis to constrain the amount of both small 
and large dust grains. The outline of the paper is as follows: 
In Section 2 we continue discussing existing constraints on 
intergalactic dust and argue for a physical model of the in- 
tergalactic dust. In Section 3 we evaluate the effects of inter- 
galactic dust on QSO colours and in Section 4 we examine 
how the unresolved SXB can put further constraints on these 
dust models. We present our results in Section 5 and sum- 
marize and discuss these results and some future prospects 
in Section 6. 



2 INTERGALACTIC DUST 

Dust grains are present in our Galaxy, in the host galax- 
ies of cosmological sources, in galaxies along the line-of- 
sight and p ossibly in the imm e diate surroundings of the 
source (e.g. ISchlegel et all 1 19981 ; lOstman et al.ll2006l . 120081; 



iGoobarl I2008T ). In this paper, we will study the effect of 
intergalactic dust along the line-of-sight to cosmological 
sources. Extinction by interstellar dust in the Milky Way 
and nearby galaxies, such as the Small and Large Mag- 
ellanic Clouds (SMC and LMC), is usually parameterized 
using mean extinction laws, A\ = /(A, Rv) where A\ is 
the extinction for wavelength A and the reddening param- 
eter, Rv = Av/(Ab — Av), is the slope of the extinc- 
tion curve in th e optical region (e.g. ICardelli et al.l 1 19891 : 
iFitzpatrickll 19991 ). Within the Milky Way, the reddening pa- 
rameter ranges between 2 <J Rv <,6 for different sightlines 
and is usually approximated as Rv = 3.1. It is not known 
if these parametrizations are valid for an intergalactic dust 
population where the dust properties might be different from 
those in the interstellar medium. Small grains preferentially 
scatter light with short wavelengths, producing a steep ex- 
tinction law with small values for Rv- Very large grains 
would produce wavelength-independent grey extinction with 
Rv oo. 

To build up extinction curves from a generic population 
of dust grains one needs to know the composition, size dis- 
tributi on, scattering and a bsorption properties of the dust 
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19771) . hereafter MRN, found that in- 
terstellar extinction in the Milky Way is well fitted using sep- 
arate populations of bare silicate and graphite grains with a 
power-law distribution of sizes where dn/da oc a -3 ' 5 . Within 
the Milky Way, graphite grains typically range in size from 
0.005 to 1 /im while silicate grain sizes range from 0.025 to 
0.25 /urn. 

We will study dust models that have either a single 
grain size or truncated MRN size distributions, varying the 
size of the smallest and largest grains a m i n and a max . The 
intergalactic dust is assumed to be distributed homoge- 
neously with a constant comoving number density of dust 
grains, n(z). It is reasonable to assume that the dust den- 
sity traces the stellar mass density in the universe, and 
we also investigate cases where the comoving dust density 
is proportional to the integrated star formation rate den- 
sity, n(z) <x J z °° dz' ( 1+z ?)* £ ( z /) ; where we use the analyti- 



cal expression for the star formation rate, p* 
iHernauist fc Springell (|2003l ). 



derived by 



2.1 Optical depth to scattering and absorption 

The total optical depth to scattering and/or absorption by 
dust grains between redshift and z cm can be calculated as 



T (Aobs , ^cm) — 



Ho 
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where n(z') is the comoving number density of dust grains at 
redshift z' . Each dust grain is assumed to be spherical, with 
radius a and geometrical cross section a(X') = 7ra 2 Q(A') 
where the efficiencies for scattering and absorption are Q SC at 
and Qabs, the efficiency for extinction is Qext = Qscat + Q a bs, 
and A' = A b s /(1 + z). The dimensionless expansion fac- 
tor £(z') — wTIm(1 + z) 3 + £7a- For a distribution of grain 
sizes, Eq. Q] becomes 
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Figure 1. Upper panel: Scattering efficiency as function of wave- 
length, Qsca(A) for spherical silicate dust grains of radii, a = 
0.01 — 5.0/^m. Solid black lines show Q B ca = rr Bca /7ra 2 , where 
Cscat is the scattering cross section, from lLaor fc Drain j il993T ) 
and the dashed red lines show the approximation made in Eq. 
[3] Lower panel: Extinction efficiency as function of wavelength, 
Qext(A) = Qaca + Qabs f° r spherical silicate dust grains of radii, 
a = 0.01 - B.Ofim. 



where j^da denotes the comoving number density of dust 
grains with radii in the range a ± da/2. 

In this paper we focus on dust composed of silicate 
(Pgrain « 3.2 g/cm 3 ) and graphite grain s (p gr ain ~ 2.3 
g/cm 3 ) (described in | Draine fe Led 11984 E aor fc Drainei 
Il993l : IWeingartner fc Draine 200 llV The origina l astronomi- 
cal silicate model constructed bv lDraine fc Led |l984) used 
laboratory measurements of crystalline olivine in the vac- 
uum ultraviolet, resulting in an absorption feature at A ~ 
0.15 /xm not seen in astronomical objects. In this paper, 
we use the smoothed astronomical silicate model dielec- 
tric function obtained by rem oving this absorption feature 
l|Weingartner fc Draind 12001). The publicly available scat- 
tering and absorption efficiencies Q sca and Q a bs are given 
for wavelengths A = 10~ 3 to 10 3 /im for grains with radii 



a = 10" 3 to 10 /xm, see Fig. [1] For wavelengths shorter than 
A <^10 -3 fim (corresponding to X-ray e nergies E ;>0.5 keV) 
we use the scattering approximations (lAlcock fc Hatchettl 



wc use the scattering approximations (Alcock & . 
1 19781 : iMiralda-Escudil 19991 : iDiikstra fc Loebll2009T l 



0.7(a/^m) 2 (6 keV/E) 2 
1.5 



Qscat < 1 

otherwise. 



(3) 

In Fig. [2] we show the corresponding total-to-selective 
extinction ratio Ry = Av/{Ab — Ay) for the graphite 
(black lines) and silicate dust (red lines) models employed 
in this paper. Solid lines correspond to single size models, 
dashed and dash-dotted lines to truncated MRN grain size 
distributions where aw n and a max are varied, respectively. 
The inclusion of large grains gives Ab ~ Av, which drives 
Ry 1 — > 0. Note that the Rv-value only refers to the mod- 
elled optical extinction in the B- and V-band and does not 
necessarily describe the entire extinction curve as observed 
for Milky Way dust. In fact, the variations in the Q-values 
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Figure 2. Ry 1 = (Ab — Ay) /Ay reddening for dust models with 
graphite (black lines) and silicate (red lines) grains with a single 
grain size (solid lines) or models with truncated MRN grain size 
distributions with 0.05/im < a m i a < l.Oyum and a max = 2.0 ^m 
(dashed lines) or a m i n = 0.02 fim (graphite) or a m i n = 0.05 fim 
(silicate) and 0.1 < a max < 2.0 fim (dotted-dashed lines). The 
inclusion of grains with radii larger than a ^0.1 — 0.2/xm gives gray 
extinction with Ry 1 — > 0. The horizontal dotted line shows the 
average reddening for interstellar dust in the Milky Way, Ry = 
3.1 



for silicate dust grains a ~ 0.25 pm can also introduce a 
small blueing-effect with Ry <^0. 



2.2 X-ray scattering by dust 

Small angle scattering of X-rays by dust grains along 
the line-of-sight can produce diffuse halos around X- 
ray point sources, as observed for many galactic X-ray 
sources. Measurements of the intensity and angular ex- 
tent of such hal os provide a quantitative test of interstel- 
lar grain models dMauche fc Gorensteinlll 986; Ma this fc Led 
1 199 it ISmith fc Dweklll998l : iDraind 120031 '). In a similar way, 
intergalactic dust grains would produce diffuse halos around 
X-ray sources on cosmological distances. In common with X- 
ray halos seen around galactic point sources, the angular ex- 
tent of an intergalactic h alo is determined by the scattering 
prop e rties of the grains (|Evans et all 1 19851 : IVaugh an et al.l 
120061 : iPetric et al.ll2006l) . 

The differential cross sections in the Rayleigh-Gans ap- 
pro ximation for a spherical dust grain of radius a is given 
by l|Havakawalll"970l : iMauche fc Gorensteid[l986l ) 
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where x — (iira/X) sin(0/2) and ji(x) = (sina;)/a; 2 — 
(cosx/x) is the spherical Bessel function of the first order. 
Ae is a normalization depending on the grain composition 
and energy E of the X-rays being scattered, 
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3 g cm -3 J \ Z 
where F(E) is the atomic scattering factor (|Henke et al.l 



4 J. Johansson et al. 



19821 Il993h , Z is the atomic charge, M is the atomic mass 
number and p gr ain is the mass density of the dust grain. 
The central core of Eq.[4]is approximately gaussian, and the 
root mean square (rms) of the scatt ering angle 6 indicates 
the t ypical size of the scattered halo l|Mauche fc Gorensteinl 
Il986ft . 6Us « 62.4 (1.0 /an/a) (IkeV/E) arcsec. 

For low optical dep ths (r B ca <S 1), the h alo fluence per 
unit scattering angle is (|Vaughan et alj|2006h . 
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where N g is the dust grain column density along the line- 
of-sight and Fx is the unscattered flux of the X-ray point 

source. 

ISmith fc Dwekl (|l99Sf ) showed that the Rayleigh-Gans 
approximation and the more exact Mie theory are in close 
agreement, given that the photon energy E (in keV) of the 
X-ray being scattered is larger than the grain radius a (in 
/im). For our purposes, where we want to study the con- 
tribution of dust scattered X-ray halos around AGN to the 
0.5-2.0 keV Soft X-ray Background, we note that the aver- 
age photon energy at the site of the scattering dust grain 
is higher than the observed energy by a factor (1 + z gra in), 
which tends to alleviate the error introduced by the use of 
the Rayleigh-Gans approximation for the larger grains. 



3 CONSTRAINTS FROM QUASAR COLOURS 

QSOs have been found to be relatively homogeneous in 
terms of colours and spectral features over a large redshift 
range . The SDSS DR7 Quasar Catalogue |Schneider et al.l 
l2010t ) contains 105783 spectroscopically confirmed QSOs 
with optical magnitudes measured through five broadband 
filters (u, g, r, i, z). We only include point sources (removing 
QSOs flagged as extended) brighter than the limiting magni- 
tudes [u,g,r,i,z] = [22.3,22.6,22.7,22.4,20.5], correspond- 
ing to S/N greater than 5:1. Objects dev iating by more than 
2a from the mean colour are r ejected l|Mortsell fc Goobarl 
120031 : 1 Ostman fc Mortsell|[2005l ). 



3.1 Composite quasar spectral template 

For our purposes it is preferable to use an unabsorbed QSO 
spectrum. We use the SDSS median co mposite spectrum 
(as derived by IVanden Berk et ail 120011) spliced with the 
HST radio-quiet composite spectrum ( Telfer et al.ll2002h to 
achieve a wavelength range 302 < A < 8552 A. The con- 
tinuum is well fitted by a broken power law (F oc A QA ) 
with spectral indices a\ = —0.24 blueward of Lya (A <4200 
A), ol\ = -1.56 for 1200 < A < 4850 A and a A = -0.42 for 
A > 4850 A, see Fig. [3] Although the statistical uncertain- 
ties in the continuum fits are small, the quoted systematic 
uncertainties in the spectral indices are cr QA ~ 0.10 — 0.15. 



3.2 Simulated Quasar Colours 

The observed dust attenuated flux F obs of an object at red- 
shift z cul observed at wavelength A D b s is given by 

-Text(A obs ,z om ) 
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Figure 3. The composite QSO spectral template together with 
a broken power law, F\ oc A aA , with spectral indices ct\ = —0.24 
(dotted line, for A < 1200A), a x = -1.54 (dashed line, for 
1200 < A < 4850 A) and a x = -0.42 (dotted-dashed lines for 
A > 4850A). Shown are also the SDSS ugriz transmission curves. 



where F cm is the intrinsic flux and the optical depth to 
extinction r cxt can be calculated using Eq. [1] and [2] for 
graphite and silicate dust grains. The dimming (i.e. in- 
crease in apparent magnitude) of light emitted at redshift 
2 em observed at a wavelength A b s is then Am (A b s , z cm ) — 
if^Tcxt (Aobs, z em ). For each dust model we simulate the at- 
tenuation of the median QSO template spectrum and per- 
form synthetic photometry using the SDSS itgrriz-filter func- 
tions. For each pair of filters X and Y, we compare the sim- 
ulated colours with observed colours, 



Apr - y) z = (x - y)° bs - (x - y)T 



(8) 



For each set of filters X and Y, we calculate the mean colour 
(X-Y)° hB in redshift bins of Az = 0.05 in the redshift range 
0.5 < z < 3. At higher redshifts the flux decrement due to 
the Lya break makes the bluer bands inefficient, so for the u- 
and g-band we only include QSOs with z < 1.5 and z < 2.3 
respectively in the analysis. 

The probabilities of different dust scenarios for each 
colour X — Y are then calculated using 

N 

X \X-Y) = J2 A(X-Y) Zt V(X-Y)-lA(X-Y) Z] , (9) 
where i and 7 index the N redshift bins. The covariance 



matrix V(X 
V(X - Y) u 



Y)i,j 



IS 



t(x - r) c 2 a i + 

<r{X - r)t 2 cmp + <5^(A -Y) 2 obB . 



(10) 



) = F e ] 



(7) 



The uncertainties in the photometric calibration a(X — 
Y)c&\ JiO.02 and the typical observed colour variance in 
the QSO sample a{X - Y) 2 obs ^0.003 for redshifts 0.5 < 
z < 2.5. The major source of uncertainty in this analysis is 
cr(X — F)temp) due to the tilting of the spectral template ac- 
cording to the systematic uncertainties of the measured con- 
tinuum spectral indices, <t Qa « 0.10. In practice, we combine 
the results from all possible colour combinations taking into 
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account the correlations between different colours and red- 
shifts. This is accomplished using Monte Carlo simulations 
to create a total covariance matrix. 



4 CONSTRAINTS FROM THE SOFT X-RAY 
BACKGROUND 

4.1 The Cosmic X-ray Background 

The cosmic X-ray background comprises the integrated 
emission of X-ray sources, primarily AGN, extended emis- 
sion fro m galaxy clus t ers, f aint starburst and "normal" 
galaxies. iMoretti et al.l (|2003h measure the total Soft X-ray 
Background (SXB) in the 0.5 - 2 keV band to be (7.53 ± 
0.35) x 10" 12 ergs/s/cm 2 /deg 2 and find that 94± 7 6 7 % of the 
SXB can be ascribed to discrete sourc e emission, i.e ~ 6 ±6% 
of the SXB is unresolved. Similarly, iHickox fe Markevitchl 
l|2007f) find that the total SXB in the 1-2 keV band is 
(4.6 ± 0.3) x 10~ 12 ergs/s/cm 2 /deg 2 . After excluding re- 
gions of radius <^2 — 20" around detected X-ray, HST and 
Spitzer sources, they find a remaining unresolved SXB of 
(3.4 ± 1.4) x 10~ 13 ergs/s/cm 2 /deg 2 , which is 7.3±3.0% of 
the total SXB. 



4.2 AGN X-rays scattered by dust 

AGN produce a dominant fraction (~ 80%) of the SXB at 
energies 0.5 < E < 2.0 keV. If dust pervades throughout 
the intergalactic medium, its scattering opacity would pro- 
duce diffuse X-ray halos around AGN. If the angular extent 
of these halos are large enough, they will effectively con- 
tribute to the unresolved SXB. This fact, along with the 
observational upper limits on the unresolved SXB, allows us 
to c onstrain the amount and properties of intergalactic dust. 

iDiikstra fc Loebl (|2009h argue that dust scattered X- 
ray halos around AGN can maximally account for a fraction 
/halo ~ 5 — 15% of the SXB, allowing them to put upper 
limits on the opacity of intergalactic dust grains Td us t(z = 
1,A = 8269A) jCO.l - 0.2(/ ha i o /10%), depending on grain 
size or size distribution. In the spirit of this approach, we 
calculate the flux of X-ray photons observed in the Soft X- 
ray energy range E = 0.5 — 2.0 keV, which is expected to be 
scattered into halos, 

fhaio = / S{E)dE / T{z') x e 7 **^'' > - 1 \dz', 

Jo.Skcv JO 1 J 

(11) 

where the optical depth to scattering r SC a is calculated using 
Eq. Q] and [5] We take the observed spectral energy density 
of the AGN to be given by S(E) cx E~ 1A normalized to 
unity over the energy range (the final results of the analy- 
sis depend only weakly on the choice of th e spectral index, 
IMoretti et aill2003l ; IDiikstra fc Loebll2009h and, 

™ = OTWWY (12) 

The comoving X-ray emissivity L(z) of AGN is expressed in 
terms of an integral over the AGN luminosity function 

£(z) =j n " X L x i>{L,z)d\o%L. (13) 



1.05 r 




Figure 4. Fraction of halo flux falling outside an aperture of 
radius r = rgo = 2.2" (blue lines), r = Argo (green lines), r = 
9rgo (red lines) calculated using Eq. (0 integrated over angles 
and an energy range between -E m i n = 0.5 keV and -Emax = 2.0 
keV assuming a spectral energy distribution, S(E) oc E' 1 . Solid 
lines show results for single grain sizes, dashed lines show models 
with truncated MRN grain size distributions with 0.01 < a m i n < 
1.0 (im and a ma x = 2.0 /an and dotted-dashed lines show models 
with truncated MRN grain size distributions with a m - ln = 0.05 ^m 
and 0.1 < Omax < 2.0 yum. 



The X-ray luminosity functio n if)(L, z)dlog L is de scribed 
by the fitting formula of iHopkins et aD ([20071 ) with 
logLmin = 40.4 and logL max = 48.0. The total contribution 
to the SXB in the observed 0.5-2.0 keV band is obtained 
by integrating the X-ray emissivity of AGN over redshift, 
Jagn = J °° T{z')dz = 5.74 x 10~ 12 ergs/s/cm 2 /deg 2 . 

In the following, we conservatively argue that the 
amount of flux from AGN being scattered far from the orig- 
inal line of sight, i.e. into an extended halo, Fhaio, will effec- 
tively be measured as a diffuse background and can thus not 
be larger than the measured value of the unresolved SXB, 

fuSXB- 

When IHickox fc Markevitchl (|2007f ) measure the unre- 
solved part of the SXB, they exclude regions of radius 
8 < rgo — 2.2 arcsec for HST and Spitzer sources unde- 
tected in X-rays and regions 8 < 4 — 9rgo for X-ray detected 
sources. Dust grains having a Jj2 /jm will scatter X-rays into 
very compact halos with angular sizes comparable to these 
excluded regions, thus contributing less flux to the unre- 
solved SXB. While the scattering is dominated by the largest 
grains, inclusion of smaller grains would cause the halo flux 
to fall off less rapidly at large angles, and the resulted tail 
could contain a substantial flux which would be effectively 
be considered as an unresolved background. We verify that 
a major fraction, g(8 > r), of the dust scattered halo flux 
falls outside these apertures (see Fig. [4} for dust models with 
single grain sizes or MRN grain size distributions, and will 
henceforth assume that si milar exclusion regio ns apply to 
the SXB measurements of IMoretti et al.l (|2003h . 

We calculate the probability of different dust scenarios 
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using 



XSXB — 



/ pobs 
(,-^uSXB 



mod\ * 
halo / 



(14) 



where F^sxb is the unresolved Soft X-ray Background in 
the 0.5-2 keV band as calibrated by iMoretti et all (|2003f ) 
and F£1°q — F ha i • g(6 > r) is the modelled halo flux falling 
outside an aperture of radius r. For values i*2a ^ ^uSxb, 
we set the corresponding x 2 -value to zero. The model un- 
certainty will be negligible compared to the uncertainty of 
the unresolved SXB a obs = 4.52 x 10~ 13 ergs/s/cm 2 /deg 2 . 



5 RESULTS 

We present our results from the combined analysis of QSO 
colours and the unresolved SXB together with the expected 
dimming in the restframe B- and F-band for a source at 
redshift z = 1, Ab(z = 1) and Ay(z = 1). We also express 
our constraints as upper limits on the density parameter, 
f^dust = (n^a 3 p Eraill ) /pent, where p clit = 1.88 x 10~ 29 ft 2 g 



'dust — [n^-a 3 Pgr, 

cm" 3 . 



5.1 Dust models with single grain sizes 

In Fig. O constraints from the QSO colour analysis for the 
single size graphite dust model are shown. Regions (from 
yellow to red) indicate allowed dust models at 68%, 90%, 
95% and 99% confidence levels. Black lines show the rest- 
frame B- and V-band extinction in magnitudes for a source 
at z — 1, Ab(z = 1) and Av(z = 1). Small size dust grains 
are more effectively constrained and the QSO colour analy- 
sis rules out any dimming Ab(z = 1) J>0.05 mag for grains 
with a <J0.2 pm. As is evident from Fig. [5l grains with radii 
a ;>0.3 pm produce very little reddening for cosmological 
sources and Ab{z = 1) ~ Ay(z = 1) ,>0.4 mag can not be 
excluded. 

Complementary to the QSO colour analysis, the SXB 
analysis (see Fig. [S| is most effective for grains with sizes 
a ;>0.2 pm, since for these grain sizes the scattering cross 
section for X-rays <7 sca tx a 4 . Here, the unresolved SXB mea- 
surements constrain the the restframe B-band extinction, 
Ab(z = 1) <J0.25 mag at the 99% confidence level. Com- 
bined constraints are shown in Fig. [7] For a dust population 
of single size silicate grains, results are both qualitatively 
and quantitatively similar as can be seen from the QSO and 
SXB combined constraints shown in Fig. [8] 

Assuming the IGM consists solely of dust grains smaller 
than a ~ 0.2 pm with a constant comoving number density, 
we can constrain the dust density parameter fidust <J10~ 6 — 
10- 5 (p grai n/3 gem ), which is close to the astrophysically 
interesting levels implied by studies of stellar evolution and 
metallicity as a function of redshift. For grains larger than 
a ~ 0.2 pm, the constraints on fidust are degenerate with a. 



5.2 Dust models with continuous grain size 
distributions 

Any realistic dust model will have a distribution of grain 
sizes. Here, we use an MRN model of separate populations 
of bare silicate and graphite grains with a power-law distri- 
bution of sizes where dn/da oc a -3 ' 5 . 



First, we study the impact of the smallest grains in 
the MRN size distributions by varying o m in and keeping 
flmax = 2.0 pm. Figures [TT] and rj2] show the combined 
QSO and SXB constraints on silicate and graphite grains 
with MRN size distributions and a constant comoving num- 
ber density. These dust models all have a reddening pa- 
rameter 6 <,Rv <J°°, representing very grey dust. Allow- 
ing dust grains a <J0.1 pm, produces enough reddening at 
optical wavelengths that the QSO colours can limit the 
extinction in the B-band, Ab(z — 1) <J0.10. For mod- 
els with a m i n !>0.1 — 0.3 pm (completely grey dust) the 
SXB data helps to limit the extinction in the B-band to 
A B (z = 1) ^0.25. 

An interesting aspect of the MRN distribution is that 
the opacity is dominated by grains with small radii whereas 
the total dust mass is dominated by the grains with large 
radii. Thus, removing the very small grains can affect the 
opacity dramatically, without radically changing the total 
mass in dust. Therefore, for MRN distributions with a max = 
2.0 pm, the combined upper limits on the dust density are 
quite large, fi dust <;10~ 5 - 10 _4 (p gra i n /3 g cm -3 ). 

To relate to the detection of dust reddeni ng of QSOs out 
to larg e distances around z ~ 0.3 galaxies by I Menard et al] 
(2010), we focus on dust models with a reddening parameter 
close to their measured value of the reddening parameter, 
Rv = 3.9 ±2.6. Assuming a constant comoving dust density 
and extrapolating their result to z — 1 yields a lower limit on 
the extinction in the restframe B-band of Ab(z = 1) ^0.02 
magnitudes. If we keep a m i n fixed at 0.05 pm (silicate grains, 
Fig. O or at 0.02/^m (graphite grains, Fig. HOjl and vary the 
largest grain sizes 0.1 < a max < 2 pm, the dust models span 
a range of Rv between ~ 1.3 — 7. For these dust models, 
the combined QSO colour and SXB analysis puts an upper 
limit on the extinction in the restframe B-band, Ab(z = 
1) <,0.05 — 0.10 mag and correspondingly limits the dust 
density parameter fWst <T0~ 6 — 10~ 5 (p gra i n /3gcm -3 ). 



6 SUMMARY AND DISCUSSION 

In this paper, we studied the optical extinction and X- 
ray scattering effects of intergalactic dust grains. We find 
that dust distributions including graphite and silcate grains 
smaller than a <C0.1 — 0.3 pm, produce enough reddening 
(1.3 <Rv <10) to rule out extinction in the restframe B- 
band, Ab(z — 1) J>0.05 — 0.10 mag. By combining the con- 
straints from the QSO colour analysis with constraints from 
the unresolved SXB, we are able to rule out any systematic 
dimming in the restframe B-band for a source at redshift 
z = 1 of Ab(z — 1) ;>0.25 mag for a wide range of grain 
sizes (and consequently a large range of reddening param- 
eters, 1.3 <,Rv ^°°)- These results have been derived as- 
suming a constant comoving dust density. In the case that 
the comoving dust density is proportional to the integrated 

the 



star formation rate density, n(z) oc J z °° dz' ( 1+z ?)* £( y ) , 
dust density decreases with redshift, making the universe 
increasingly transparent relative to the constant comoving 
dust density models. Results from these models are quali- 
tatively similar to the models with constant comoving dust 
densities, except that the constraints on the induced dim- 
ming for sources at redshift z = 1 (see Fig. ll3[) are somewhat 
weaker, excluding Ab(z = 1) !>0.05 for small grains with 
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a <0.1 /im and Ab{z = 1) >0.35 mag for large grains with 
a J>0.2 /im. However, for the same reasons, the upper lim- 
its on the dimming of sources at higher redshift (z J>2) are 
stronger compared to the models with constant comoving 
dust density. 

Since the error budget is dominated by systematic tem- 
plate uncertainties, the increas e in the number of QSO s 
doe s not improve the limits fro m I Mortsell fc Goobarl l|2003l ) 
and lOstman fc Mortselll (2005), although the results are not 
directly comparable since the dust models probed are not 
equivalent. However, we regard the current error treatment 
as more realistic and also more effective in handling all pos- 
sible colour combinations while accurately taking into ac- 
count all possible colour and redshift correlations. In order 
to improve the current results from QSO colours (which are 
most effective in constraining the density of small size dust 
grains) it would be necessary to be able to decrease the sys- 
tematic uncertainty of the continuum slope of the template 
spectrum. For larger dust grains, where the reddening ef- 
fects are smaller and the dimming and density constraints 
from the unresolved SXB are weaker, limits could be im- 
proved by a better understanding of the different contribu- 
tions to the unresolved SXB. Furthermore, improved con- 
straints may be obtained by stacking observations of X-ray 
point sources. The absence of a single X-ray h alo around a 
z = 4 .3 QSO observed with Chandra, allowed IPetric et al.l 
(2006) to place an upper limit on fidust <^2 x 10 -6 assuming 
a constant comoving number density of dust grains with a 
characteristic grain size a ~ 1/im. This result is relaxed in 
ICorrales fc Paerels! (|2012h . who find that Q. dust ^10~ 5 for a 
wider range of dust models. However, an increased number 
of sources will allow us either to measure or to improve the 
constraints on large dust grain densities, for which the cor- 
responding limits on the dimming of cosmological sources 
are currently the weakest. 

It should be noted that while smaller grains 
(a <;0.05 fim) may be either destroyed by sputtering or 
unable to travel far from formation sites as they are in- 
efficiently pushed away by radiation pressure, large grains 
(a J>0.25 ^m) may be too heavy and remain trappe d in the 
gravi t ational field of the g a laxy where they formed jAguirrel 
1999:; lAguirre et all l200ll : iBianchi fc Ferraral 120051). Com- 
bined with stellar evolution models (lFukugitall201ll ). the 
lower limits obtained in iMenard et all "! 2010l ) and the up- 
per limits derived in this paper, we conclude that an MRN 
dust distribution of either silicate or graphite grains - or 
possibly an admixture thereof - with a m i n ~ 0.05 /im and 
~ 0.25 [iva and f^dust ~ & few xlO 6 (p g rain/3gcm 3 ), 
constitute a viable intergalactic dust model. 
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QSO constraints 
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Figure 5. Constraints from the QSO colour analysis on dust 
models with single size graphite grains and a constant comoving 
number density. Regions (from yellow to red) indicate allowed 
dust models at 68%, 90%, 95% and 99% confidence levels. Black 
lines show the restframe B- and V-band extinction in magnitudes 
for a source at z = 1, Ag(z = 1) and Ay(z = 1). 



SXB constraints 
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Figure 6. Constraints from the unresolved Soft X-ray Back- 
ground on dust models with single size graphite grains and a 
constant comoving number density. Regions (from yellow to red) 
indicate allowed dust models at 68%, 90%, 95% and 99% confi- 
dence levels. Black lines show the restframe B- and l^-band ex- 
tinction in magnitudes for a source at z = 1, Ag(z = 1) and 
A v (z = 1). 



Vaughan S., Willingale R., Romano P., Osborne J. P., Goad 
M. R., Beardmore A. P., Burrows D. N., Campana S., 
Chincarini G., Covino S., Moretti A., O'Brien P. T., Page 
K. L., Supper M. A., Tagliaferri G., 2006, ApJ, 639, 323 

Weingartner J. C, Draine B. T., 2001, ApJ, 548, 296 



Combined QSO and SXB constraints 
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Figure 7. Combined constraints from the QSO colour analysis 
and the unresolved Soft X-ray Background on dust models with 
single size silicate grains and a constant comoving number den- 
sity. Regions (from yellow to red) indicate allowed dust models at 
68%, 90%, 95% and 99% confidence levels. Black lines show the 
restframe B- and V-band extinction in magnitudes for a source 
at z = 1, Ag(z = 1) and Ay(z-l). 
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Figure 8. Combined constraints from the QSO colour analysis 
and the unresolved Soft X-ray Background on dust models with 
single size silicate grains and a comoving number density propor- 
tional to the integrated star formation rate. Regions (from yellow 
to red) indicate allowed dust models at 68%, 90%, 95% and 99% 
confidence levels. Black lines show the restframe B- and V-band 
extinction in magnitudes for a source at z = 1, As [z = 1) and 
A v (z = l). 
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Figure 10. Combined constraints from QSO colour analysis 
and the unresolved Soft X-ray Background on dust models with 
graphite grains with a truncated MRN distribution of sizes 
a min = 0.02 /im and 0.1 < a max < 2.0/xm and a constant co- 
moving number density. Regions (from yellow to red) indicate 
allowed dust models at 68%, 90%, 95% and 99% confidence lev- 
els. Black lines show the restframe B- and V-band extinction in 
magnitudes for a source at z = 1, Ag(z = 1) and Ay(z = 1). 
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Figure 9. Combined constraints from QSO colour analysis and 
the unresolved Soft X-ray Background on dust models with sil- 
icate grains with a truncated MRN distribution of sizes a m ; n = 
0.05 (im and 0.1 < a max < 2.0 /im and a constant comoving num- 
ber density. Regions (from red to yellow) indicate allowed dust 
models at 68%, 90%, 95% and 99% confidence levels. Black lines 
show the restframe B- and V-band extinction in magnitudes for 
a source at z = 1, Ab(z = 1) and Ay(z = 1). 
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Figure 11. Combined constraints from the QSO colour anal- 
ysis and the unresolved Soft X-ray Background on dust mod- 
els with silicate grains with a truncated MRN distribution of 
0.01 < Omin < 1 M m an d "mm = 2 fim and a constant comov- 
ing number density. Regions (from yellow to red) indicate allowed 
dust models at 68%, 90%, 95% and 99% confidence levels. Black 
lines show the restframe B- and V-band extinction in magnitudes 
for a source at z = 1, Ab(z = 1) and Ay(z = 1). 
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Figure 12. Combined constraints from the QSO colour analysis 
and the unresolved Soft X-ray Background on dust models with 
graphite grains with a truncated MRN distribution of 0.01 < 
a min < 1 A 1111 an d «max = 2 /im and a constant comoving number 
density. Regions (from yellow to red) indicate allowed dust models 
at 68%, 90%, 95% and 99% confidence levels. Black lines show the 
restframe B- and V-band extinction in magnitudes for a source 
at z = 1, Ag(z = 1) and Ay(z = 1). 
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Figure 13. Combined constraints from the QSO colour analysis 
and the unresolved Soft X-ray Background on dust models with 
single size graphite grains and a comoving number density propor- 
tional to the integrated star formation rate. Regions (from yellow 
to red) indicate allowed dust models at 68%, 90%, 95% and 99% 
confidence levels. Black lines show the restframe B- and V-band 
extinction in magnitudes for a source at z = 1, Ag(z = 1) and 
A v (z = l). 



